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Abstract 

This paper proposes a unified approach that enables the Wishart distribution to 
■ be studied simultaneously in the real, complex, quaternion and octonion cases. In 

ON ' particular, the noncentral generalised Wishart distribution, the joint density of the 

eigenvalues and the distribution of the maximum eigenvalue are obtained for real 
normed division algebras. 

2 ■ 1 Introduction 

O 

Many results first described in statistical theory are then found in real cases, and the version 
for complex cases is su bsequently studi ed. This has b e en de s cribed in vari ous papers, see 
Br avais (1846 ) (cited by IWoodinef (Il956h ) and IWoodingj (|l956l) ; IJamesl (|l964 Sections 4 and 
ldF" 



X' 8): lMuirheadl (Il982h and iRatnaraiah et all (|2005l) . among many other examples. 



Using some concepts and results derived from abstract algebra, it is possible to propose 
a unified means of addressing not only real and complex cases but also the quaternion and 
oct onion cases. Part o f this a ppro ach has been use d for some time in random matrix theory, 
Edelman and Rao (|2005l) and lForresten (|2009l ). 



see 



For the sake of completeness, in the present study the case of octonions is considered, but 
it should be noted that many results for the octonion case can only be c onjectured, because 
there remain many unresolved theoretical problems in this respect, see Drav and Manoguei 



(1999). Furthermore, the relevance of the octonion case for understanding the real world 



has yet to be clarified, see iBaea (|2002i ) . 

The rest of this paper is structured as follows: Section [5] reviews some definitions and 
notation on real normed division algebras. Some results for Jacobians and generalised hy- 
pergeometric functions, together with an extension of one of the basic properties of zonal 
polynomials (which is also valid for Jack polynomials) are also given. Section [3] then derives 
the noncentral generalised Wishart distribution, and as a corollary the noncentral inverse 
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generalised Wishart distribution is obtained. In Section [H we obtain the joint density func- 
tion of the eigenvalues and the distribution of the maximum eigenvalue, the latter under 
a matrix multivariate normal distribution. All these results are obtained for real normed 
division algebras. 

2 Preliminary results 

Let us introduce some notation and results that will be useful. 



2.1 Notation and real normed division algebras 



A detailed discussion of real normed division algebras may be found in iBaea ((2002) . For 
convenience, we shall introduce some notations, although in general we adhere to standard 
notations. 

For our purposes, a vector space is always a finite-dimensional module over the field 
of real numbers. An algebra J is a vector space that is equipped with a bilinear map 
m '■ -5 x -5 -5 termed multiplication and a nonzero element 1 e 5 termed the unit such 
that m(l, a) — m(a, 1) = 1. As usual, we abbreviate m(a, b) = ab as ab. We do not assume 

5 associative. Given an algebra, we freely think of real numbers as elements of this algebra 
via the map lj t— >■ ul. 

An algebra 5 is a division algebra if given a, b £ $ with ab = 0, then either a = or 

6 = 0. Equivalently, $ is a division algebra if the operation of left and right multiplications 
by any nonzero element is invertible. A normed division algebra is an algebra J that is 
also a normed vector space with ||a6|| = ||a||||fr||. This implies that 5? is a division algebra 
and that 1 1 1 1 = 1 . 

There are exactly four normed division a lgebra s : real numbers (3?), complex numbers 



(£), quaternions ($)) and octonions (D), see iBaed (|2002l ). We take into account that 5ft 



£, .f) and O are the only normed division algebras; moreover, they are the only alternative 
division algebras, and all diy ision algebras have a real dimension of 1,2,4 or 8, which is 
denoted by ft, see iBaezI ( 2002 , Theorems 1, 2 and 3). In o ther branches of mathematics, the 



parameter a = 2/ ft is used, see lEdelman and Raol (l2005h 



Let £^ „ be the linear space of all n x m matrices of rank to < n over $ with to distinct 
positive singular values, where $ denotes a real finite- dimensional normed division algebra. 
Let 5™ xm be the set of all n x to matrices over J. The dimension of $ nxm over 5ft is j3mn. 

T 

Let A G £ nxm ; th en A* = A denotes the usual conjugate transpose. 

The set of matrices Hi e y» xm suc h that H*Hi = I m is a manifold denoted V^ n , is 
termed the Stiefel manifold (Hi is also known as semi-orthogonal (ft = 1), semi-unitary (0 = 
2), sem i-symplectic ((3 = 4) and semi- exceptional type (fi = 8) matrices, ^ee l Prav a.nd Mano.u 3 
(|l999f )1. The dimension of „ over 5ft is [f3mn — m(m — l)/3/2 — m]. In particular, m 



with dimension over 5ft, [m(m + l)/3/2 — to], is the maximal compact subgroup il^(m) of 
jC^ m and consists of all matrices H £ J mxm such that H*H = I m . Therefore, il^(m) is the 
real orthogonal group 0(m) (ft = 1), the unitary group U(m) (j3 = 2), compact symplectic 
group Sp(m) (ft — 4) or exceptional type matrices Oo(m) (ft — 8), for J = 5ft, £, fj or D, 
respectively. 

We denote by the real vector space of all S <E g" mxm such that S = S*. Let 
^P^j be the cone of positive definite matrices S £ y" xm ; then is an open subset of 
6^. Over 5ft, consist of symmetric matrices; over £, Hermitian matrices; over $), 
quaternionic Hermitian matrices (also termed self-dual matrices) and over O, octonionic 
Hermitian matrices. Generically, the elements of 6^ are termed Hermitian matrices, 
irrespective of the nature of The dimension of over 5ft is [m(m — l)ft + 2] /2. 
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Let be the diagonal subgroup of C^ nm consisting of all D G 5 raxm , D = diag(di, . . . , d m ). 

For any matrix X G yx"^ <^x denotes the matrix of differentials (dxij). Fi nally, we 
define the measure or volume element (dX) when X G $ mx ™ ; 35^ or n; see lDimitrTu 
(|2002h . 

If X G J" xm then (dX) (the Lebesgue measure in g" nxm ) denotes the exterior product 
of the /3mn functionally independent variables 



P 

(fc) 

8=1 3=1 k=l 



(dX) = A A dxij where dxy = A d:z^ 



(dS) 



If S G 6^ (or S G 2^ (to)) then (dS) (the Lebes gue measure in or in *X^(to)) denotes 
the exterior product of the m(m+ l)/3/2 functionally independent variables (or denotes the 
exterior product of the m(m — l)/3/2 + n functionally independent variables, if sn € -ft for 
all i = 1, . . . , to) 

m P 

AA*S } > 

m m (3 

f\ds u AA d 4 fe) - ifs "e3fi. 

i= 1 i<j fc=l 

The context generally establishes the conditions on the elements of S, that is, if Sij E 5ft, 
€ £, € or € 0. It is considered that 

i<_7' fc=l 2—1 i<j fc— 1 

Observe, too, that for the Lebesgue measure (dS) defined thus, it is required that S € ^3^, 
that is, S must be a non singular Hermitian matrix (Hermitian definite positive matrix). 

If A G 2)^ then (dA) (the Legesgue measure in 33^) denotes the exterior product of the 
(3m functionally independent variables 

n /3 

(fc) 



(dA) = f\ f\ dX 



i=l k=l 



If Hi G V£ „ then 



(H*dHi) = /\ A h J dhi " 

i=l j—i+l 

where H = (Hi|H 2 ) = (hi, . . . , h m |h m+ i, . . . , h„) G il^(m). It can be proved that this 
differential form does not depend on the choice of the H2 matrix. When to = 1; vf n defines 
the unit sphere in $ n . This is, of course, an (n — l)f3- dimensional surface in $ n . When 
m = n and denoting Hi by H, (H*dH) is termed the Haar measure on if 3 (to). 
The surface area or volume of the Stiefel manifold n is 

Vol (<J = / (HtdHi) = (1) 

•/H!6v£.„ T, n [np/2\ 

where r^Ja] denotes the multivariate Gamma function for the space 6„, and is defined by 
I* [a] = / etrj-AllA^-^-^^^-^dA) 



m 

m("»-iW4TJr[a-(i-l)j9/2], 



7T 

i=l 
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where etr(-) = exp(tr(-)) , | • | denotes the determinant and Re(a) > (m — l)/3/2, see 
Gross and Richards! |l987t) . 



2.2 Jacobians and and some results on integration 

Fi rst, we summar ise diverse Jacobians in terms of the /? parameter, some based on the work 
of Di mitriu (l2002h. wh ile other results are proposed as ext e nsions of real, c omplex or quater- 



nion cases, see 



Jamesl (Il954h. Ijamesl (|l964T ). iKhatril (|l965| ). iMethal (|l99l| ). rRatnaraiah et al. 



(|2005h and lLi and Xud (120091) 

Lemma 2.1 (Singular value decomposition, SVZ)). Let X 6 £^ n , such that X = ViDW* 
with Vi G V£ „, W G U^(m) and D = diag(di, ■ • • , d m ) G 2)^, di > • • • > d m > 0. J7ien 

m m 

(dX) =2" m 7r T JJdf (rt " ro+1)_1 JJ(d? -^(dD)(VIdVi)(WdW), (2) 

i — 1 i<j 

0, 0=1 

-m, /3 = 2 
-2m, /3 = 4 
-4m, /3 = 8. 

Lemma 2.2 ( Spectral decomposition). Let S £ *# m - Then the spectral decomposition can 
be written as S = WAW, where W G ii^(m) and A = diag(Ai, . . . , A m ) G £>* n , with 
Ai > ■ ■ • > A m > 0. Then 



(dS) = 2- ,n n T JJ(A, - Xjf(dA)(W*dW), 



where t is defined in Lemma \2.1\ 

Lemma 2.3. Let X G C m>n , and S = X*X G q}^. T/ien 

(dX) = 2- m |S|^ ( "- m+1)/2 - 1 (dS)(V*dVi). 
Theorem 2.1. Le£ S G TAen ignoring the sign, if Y = S 

(dY) = |s|-/9C"*— i)-3(ds). 
Consider the following property of Jack polynomials. 
Lemma 2.4. Z/X G £^ n , #ien 



(3) 



(4) 
(5) 



(6) 



where C^(B) are i/ie Jacfc polynomials of weight k of B G 6^ correspon ding to t he par - 
tition k = (hi, . . .km) of k, k\ > • • • > k m > w/it/i J^Li fcj = fc, see Iffawt/eTl 1(1 99$ ). 
Gross and Richards] 1(1 98$ ): and [a]^ denotes the generalised Pochhammer symbol of weight 



K, defined as 



[af K = f[(a-{i-l)j3/2) kl , 



where 5J(a) > (to — l)/3/2 — fe m and (a), = a(a + 1) • ■ ■ (a + i — 1). 
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Proof. See 



Diaz-Garcia and Gutierrez- Jaime j ( 2010f ). □ □ 



Now, we utilise the complexification = 6^ + i&m of 6^. That is, consist 

of all matrices X £ (J E ) raXm of the form Z = X + iY, with X, Y £ 6^. We refer 
to X — Re(Z) and Y = Im(Z) as the real and imaginary parts of Z, respectively. The 
generalised right half-plane <1> = ffiff, + i&J L in 6^ consists of all Z g <S £,' £ such that 
Re(Z) £ see lGross and Richards! (|l987l p. 801). Also, as in lDavid (|l980h . consider the 
following notation, 

oo oo oo 

E E -EEEEE- 

fc,i=0 k,5;4>£k-8 k=0 1=0 K 8 <j>£n-8 



cjf^' ^A. B) denotes th e invariant poly nomials, which are defined in Chikuse and Davisl 



()l979h and|Davisj (|l980h in the real case. iDiaz-Garcial (|2009l ) studied these invariant poly- 
nomials and many of their basic properties for real normed division algebras. 

Theorem 2.2. Let A £ * then 

I | Xr -(— D/3/2-1 etr{-XA} p if (01, . . . , a p ; h, . . . , b q ; BX)(dX) 

Jo<X<A 



r^[a]rf n [(m-l)/3/2 + l] 
r&[a+(m- 1)0/2+1] 



^.....[ajg tfj«- 4 C7j««.«(-AA l BA) 



\5>jS fcm[6^,...,[a^ 



+ (m- 1)0/2 + 1] 



where 9f ' f "' <5 is defined in Diaz- Garcid \200A . eq. (52)), see also lDaviA kl98(\). Also, „Fn de 



notes the h yper geometric function defi ned in terms of Jack polynomials, see Gross and Richards] 
198j ) and \Koev and EdelmaA kOOd ). 



Proof. This follow s immediately, expanding p F@ in terms of Jack polynomials and using 



Diaz-Garrial (|2009l eq. (5.38)) 



□ 



3 Wishart distribution 



Recall that X £ £^ n has a matrix multivariate elliptically contoured distribution for real 
nor med division algebras if its density, with respect to the Lebesgue measure, is given by 



(sec 



Diaz-Garcia and Gutierrez- Jaimezj ( 20091) ): 



/x(X) 



|£|/3n/2|0|/3m/2 



h{tr [S- 1 (X-/^)*0" 1 (X-/x)]} 



where fi £ £^ „, £ £ tyf^, £ The function h : $ — > [0, oo) is termed the generator 
function, and it is such that j^p u l3nm ~ 1 h(u 2 )du < oo and 



T[f3mn/2} 
C P (m,n) = 27rfjmn/2 




(inm — l 



h(u 2 )du 



Such a distr ibution is denoted by X ~ £f (/j, S, 0, h), f o r the real case see Fang and Zhand 
(|l990h and iGupta. and Vargal jl993l ) and iMicheas et all (2006) for the complex case. Ob- 
serve that this class of matrix multivariate distributions includes normal, contaminated 
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normal, Pearson type II and VI, Kotz, Jensen-Logistic, power exponential and Bessel distri- 
butions, among others; these distributions have tails that are more or less weighted, and/or 
present a greater or smaller degree of kurtosis than the normal distribution. 

Theorem 3.1. Let S = X*0 _1 X £ S is said to have a generalised Wishart distri- 

bution for a real normed division algebra, this fact being denoted as S ~ QW^fn, X, O, h). 
Moreover, its density function is 

WK") | 3 |«^i)M V k(2k) (tr^S + ^fo^S) (?) 

where Q = X _1 /i*0 _1 /i and h^(-) is the jth derivative of h with respect to v = trS _1 S. 
Proof. Let S = X*0 *X = Y*Y, where 

Y = ©-^X ~ ^ Xro (0- 1 /2 At; £, i mi h ) t 

with (0 1 / 2 ) 2 = 0, and so 

MY) = ^#MtrS- 1 (Y-^r(Y- /i) ] 

Let us now consider the singular value decomposition of matrix Y = ViDW*. Then by 
Lemma 12. 1[ the joint density function of Vi , D and W is 

rn rn 
z—1 i<cj 



|e|/W2 

x h [tr (£- x WD 2 W* +n)+ tr (-2/j y S- 1 WDV;)] 

x (W*dW)(dD)(V3Wi), (8) 

where /Zy = -1 / 2 /^ and £ _1 /i*0 _1 /i. Let us now assume that h can be expanded in series 
of power, that is 

fc=0 

Hence, considering only /i in ([5J 

^ M fe ) [tr (E- 1 WD 2 W* + «)1 , , , 
= E — « ( tr (-2M Y S- 1 WDVt)) fc . 

fe=0 

And from Lemma |2~^1 noting that ^ is zero for all odd k, 



L 



[tr (-2tr/x Y S- 1 WDVJ)] 2fe (VJdVi) 



omJmn/2 f I "\ 4*: 
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Observing that (|) . 4 fe / (2fc)! = 1/fc!, the joint density function of D and W is 



7r ^m„/2+r C -/3 (m) „j TJ ^(n-m+l)-l jj ((f 2 _ ^0 
z— 1 i<j 

r m [/3n/2] |S|^/2 

~ [tr (s- 1 wd 2 w* + n)] ^ c K (ns- x wD 2 w*) 



X 



fc=0 

x (W*dW)(dD). 



\J3n/2]£ 



and 



Finally, let S = Y*Y = WD 2 W*. The desired result is obtained from lemmas 
noting that (dD) = 2-" i |D 2 |- 1 / 2 (dD 2 ) and = l s l- □ □ 

Distribution ([7]) was found by iDi'az-Garci'a and Gutierrez- Jaimea 1 ( 2009 Hot the real case 
and for the general central case by Diaz-Garcia and Gutierrez-Jaimea ( 20091 ) . 

Corollary 3.1. Assume that X is a matrix multivariate normal distribution for real normed 
division algebras. Then S = X*@ X X has a Wishart distribution for real normed division 
algebras and is denoted as S ~ W^(n, E, fi). Moreover its density is 



(2/^)^«/2rP t [/3n/2]|E|/ 3 »/ ; 



| S |0(n-m+l)/a-l etT {_p ( S -1 S + n ) / 2 } 

x if (/V2;/3 2 nE- x s/4). (9) 



Proof. This follows from (J7J taking into account that for the normal case, h(u) = exp{— /3u/2} 
and C !i (m,n) = (2vr//3)-' 3 " m / 2 . □ □ 

Th i s result has be en found by Hcrz (1955) and lJamesI (|1961I ) for the real cas e: bvlJamea 
1964 ). Khatri (1965) and Ratnaraiah et al. ( 2005 ) for the complex case and by lLi and Xue 



20091 ) for the central quaternion case, among other aut hors. The general central case of ^ 



was found bv lDiaz-Garcia and Gutierrez- Jaimezl ( 2009f ). 

From Theorems l2.ll and l3.ll it is straightforward to obtain the distribution of S _1 , termed 
the inverse generalised Wishart distribution. 



(10) 



Corollary 3.2. In Theorem \3.1\ we define W = S . Then its density function is 



Tr 0mn / 2 CP(m,n)\W\ ( 



n+m+l)/2-3 



r^ l [/3n/2]|E|' 3 "/2 



E 

fc=0 



h (2k) ( trS -i w -i + n ) ci (ns^w- 1 ) 



The density (| 1 0[) was found bv lip. Wong, and Liul (|2007f ) in the real case. 



4 Eigenvalue densities 

In this section we find the general joint density function of the eigenvalues of S and the 
density of A max for the normal case. 
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Theorem 4.1. Assume that S ~ 5W^,(n, S, SI, h). Then the joint density of eigenvalues 
Ai, . .., A TO > 0, o/S is 



i—1 i<j 



ri[/3n/2]Ti[pm/2]m^/2 

» _ ft(^+0(t r n)^i'" 5 rf(A)cif ]K ' 5 (s- 1 ,ns- 1 ) 

x > > 1 ^ 2: — . (11) 

Proof. Let S = WAW* the spectral decomposition of S, where W G il^(m) and A = 
diag(Ai, . . . , A m ) G S 1 (m), Ai > • • ■ > A m > 0. Then by ([7]) and Lemma r2.2[ the marginal 
density function of A is 



2 -m^m n /2+r C ^ TO] „j J| ^(n-m+l)/2-l jj^ _ ^y, 



i—1 i<j 



00 n 



r^[/3n/2]|S|/3«/2 fe![/3n/2] £ 
x [ h {2k) (trE^WAW* + Q) Cf (nS _1 WAW*) (W*dW). (12) 

By denoting the integral in (fT2]) by J and the expanding h^ 2k ^ into series of powers, we have 

h (2k+l *> (tin) 



J = E 



n 

1=0 



00 



EE 



f (trS- x WAW*)'cf (nE -1 WAW*) (W*dW) 

ft (2*+0 (trn) 



Z! 

1=0 s 



X 



/ Cf (S-'WAW) Cl (nS _1 WAW*) (W*dW) 



2 m 7r /3m 2 /2 

ri[pm/2] 

- _ h<*H> (trft)^ g (A)Cg^ (s^nST 1 ) 
x 2^ n 7^ ■ ( 13 ) 

1=0 5;<t>£K-8 " 

The last equality is obtained by applying iDiaz-Garcia (|2009l eq. (5.1)). The desired result 
then follows by substituting ([13]) in (|T2"]). □ □ 

In the real case pT|) was obtained by IDiaz-Garcia and Gutierrez- Jaime3 ( 20091) . 

Corollary 4.1. Assume that S ~ W^(n, S, f2). TTien £/ie jomf density of eigenvalues 
Ai, . . . , A m > 0, o/S is 

m m 
^m 2 /2+r T| ^(n-m+D/3-1 JJ (A . _ ^y3 etr{ _^ n /2} 



(2//3)/ 3 W2r^[/3n/2]rf n [/3m/2]|S|^/2 

C^(A)cf ]K ' 5 (-/3S- 1 /2,/? 2 ^S- 1 /4) 



F F — ^- 



s 



Proof. The proof follows from (jlll) taking into account that for the normal case, h(u) = 
exp{-^u/2} and C p (m, n) = (2 7 r//3)-' 3m ™/ 2 from where h^ 2k +l \u) = (-/3/2) 2fc+i exp{-^u/2} 
and observing that Cjf ]K ' 5 (aA, 6B) = o*6 I Cf 1M (A,B) ) see |Piaz-Gard3 (I2009L eq. (5.8)). 

□ □ 



The result in Corollarv l4 . 1 1 was ob tained by Davis J 1980) for the real case; bv lRatnaraiah et~al. 



(2005) for the complex case, and by Li and Xud ( 20091 ) for the central quaternion case. 

Theorem 4.2. Let A G 3> and consider that S ~ W^(n, S, SI) then the probability f [S < 
A] is 

rj[(m - l)/3/2 + l]|A|^"/ 2 etr{-^/2} 
(2//3)^W2rg,[/3(n + m- l)/2+ ljjsj^P 



,.,-,„,„:,.=«., fc!«D9n/2]I + m - l)/2 + 1]J 



» x . [/3n/2]£ ^ jre,0 C^ K '°(-/3S- i A/2, y 8 2 nS- 1 A/4) 

x E E 



Proo/. The proof follows from © and Theorem |2~2"1 □ □ 

Note that A max < y is equivalent to S < yl. Therefore, the distribution of A max is 
obtained by letting A = yl in Theorem 14. 2 [ and hence: 

Corollary 4.2. Assume that S ~ W^(n, X, SI) and let A max be the maximum eigenvalue of 
S. Then %\ max < y] is 

r^[(m - l)/3/2 + l]/ m "/ 2 etr{-^r2/2} 
{2/Pf mn l 2 Tt[P{n + m - l)/2 + ljjsj^V 2 

[/3n/2]^ ^"'V+'cf "''(-^A 



x 



E E 



fe^o Mtft-* fe"!D9n/2]g [0(n + m - l)/2 + 1]J 

This latter result was found by IRatnaraiah ef aD (l2005h for the complex case. 



Conclusions 

As shown in this paper, it is possible to take a general approach to the theory of distri- 
butions, for the real, complex quaternion and octonion cases simultaneously. However, any 
generalisation entails a cost, in this case the need to use some concepts, definitions and nota- 
tion from abstract algebra. Thus, any reader interested in a particular case - real, complex, 
quaternions or octonions - need simply take the particular value of /3 in order to obtain the 
results desired. 

In summary, the noncentral generalised Wishart distribution, the joint distribution of the 
eigenvalues and the maximum eigenvalue distribution are found under a unified approach 
that allows the simultaneous study of the real, complex, quaternion and octonion cases, 
generically termed distributions for real normed division algebras. 
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